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Abstract. In this paper we discuss delta shock interaction problem for a 
pressureless gas dynamics system with two different ways of approaching the 
subject. The first one is by using shadow wave solution concept. The result of 
two delta shock interactions is delta shock with non-constant speed in a general 
case. The second one is by perturbing the system with a small pressure term. 
The obtained perturbed system is strictly hyperbolic and its Riemann problem 
is solvable. We compare a limit of a numerical wave front tracking results as 
small pressure term vanishes with the shadow wave solution. 

Key words: weighted shadow waves, delta shock waves, wave front tracking, 
Riemann problem, interactions 

1. Introduction 

Consider the one-dimensional Euler gas dynamics system given by 

dtP + dr,{pu) = 
dt{pu) + dx{pv? +p{e,p)) = 0, 

where p is the density, m ~ pu la the momentum, p{e,p) — epo(p) is the scalar 
pressure, £ << 1 and po{p) — p'' /j ■ Taking e in ([1]), we obtain the pressure- 
less gas dynamics model (PGD model in the rest of the paper), also called sticky 
particles model (in fT3]) 

dtp + d^{pu) ^ 
dtipu) + d4pu^) = 0, {x,t) eRxR+. ^' 

System ([1} can be considered as a perturbation of system ^ which is weakly 
hyperbolic with a double eigenvalue Ai = A2 = u. All entropy pairs (77, q) with 
a semiconvex function i] are given by i] := pS{u), q := puS{u), where S" > 
(the entropy function r; is semi-convex with respect to the variable {p,pu)). The 
Riemann problem 



(1) 



p{x,0) = \P'-^ , u{x,Q) = \ "<^' (3) 

has a classical entropy solution consisting of two contact discontinuities connected 
with the vacuum state (p = 0) if U{) < ui: 

{{po,uo), x<uot, 
(Ojipix/t)), uot < X < uit, 
(pi,ui), X>Ult, 

where ip{y) = y. We are now turning to the case mq > ui when there is no classical 
solution to the Riemann problem ^ [3]) . 
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Throughout this paper, the foUowing constants wiU be fixed: 

1 fe 

7=1 + 2e, 0<e<-, k= and p = . 

2 V7 

2. Elementary waves of the perturbed system 
The eigenvalues of system ([T]) are 

7-1 

Ai = U - 2 , 

and the corresponding eigenvectors are 

ri = (-1, -u + K^p^Y ^ 



(4) 



(5) 



(6) 



We have chosen an orientation such that W Xi ■ ri > 0, i = 1,2, since both fields are 
genuinely nonlinear. The corresponding Riemann invariants of system ([T]) are 



s — u + — j-^(/o'^ — 1) : 1-invariant, and 
r = u — —^(/O^ ^ 1) • 2-invariant . 

The rarefaction curves through the point (po,'"o) are given by 

u — Mo = —-^{p'^ ~ Po) I < p < Po : 1-rarefaction curve, 
u — uo — — ~ Po) 1 Po ■ 2-rarefaction curve, 

while the shock curves through the point (po,uo) are given by 



and 



u — Uq 



U ~ Uq — Ki 



I P^-P6 
PoP{p - Po) 



{p — Po), p > Po : 1-shock curve. 



(7) 



(8) 



(9) 



p-'-p'o 



{p — Po), < p < Po ■ 2-shock curve. (10) 



Popip-Po) 

With the Riemann invariants, shock curves starting from the point (ro, so) are 



Si 



ro - r = Kp„ 



So - s = kPq 



{a-l){a'' - 1) 



(a - l)(aT - 1) 



V7' 



e 
-1 



(11) 



where = r{po,Uo), sq — s{po,Uo) and a = p/po > 1, and 



S2 ■■ < 



'(l-a)(l-aT) ^~a'' 
So- s = kpq[ y h V7 



(12) 
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where ro ~ r(po,Mo), so — s{po,uo) and < a = p/po < 1- The corresponding 
rarefaction curves are given by 

Ri : r>ro,s = sq, (13) 

and 

i?2 : s> So, r = tq. (14) 

It is clear that from ((TTl [T2)) we have that rg — r > sq — s holds for Si curve and 
So — s > rg — r holds for 5*2 curve, respectively. 

The Riemann problem for system ([T]) with initial data ([3]) was solved by Riemann 
[12j , and the result is summarized in the following theorem (the proof can be found 
in Courant-Friedrichs [4] and Smoller [14]). 

Theorem 2.1. [Ti Consider system fl]) with initial data {Jj. Suppose that ui — 
< -j-^{Pi + Po); equivalently Sq — ri > — H^b/l, Then there exists a unique 
solution composed of constant states {po,Uo) = {ro,So), iPm,Ujn) = (?'m,'Sm) and 
(pi,Ui) = (ri,si) separated by centered rarefaction or shock waves satisfying the 
following estimates: 

r{x, t) = r{p{x, t),u{x, t)) > min{ro, ri}, 
s{x,t) = s{p{x,t),u{x,t)) < max{so, si}. 

The amplitude of the waves is denoted by 

P ■— fm — fo '■ amplitude of an 1-wave, , , 

X := si — Sm : amplitude of a 2-wave. 

Here /?, x > for centered rarefaction waves and /3, x < for shock waves; absolute 
values Ixl are called strengths of P and x, respectively. 

We shall use that notation throughout the rest of the paper. 

3. Local Interactions Estimates 

Our first task is to obtain a sharp estimate of wave strengths with respect to e as 
much as possible. In order to do that, we shall present some assertions from [TT] 
together with modified proofs, since certain changes in estimates will be useful for 
our investigation. 

Theorem 3.1. |11| The shock curve Si starting at the point (ro,so) is given by 

So- s = gi{ro~r,po) = I /ii(a)U=ai(/3/Kpg) r < ro, (17) 

Jo 

where < g'i{P,po) < 1 and g'{iP,po) > (0- The shock curve S2 starting at the 
point (ro, sq) is 

So-s 



ro-r ^ g2{so- s,po) ^ i /i2(a)U=a2(x/'^pg) ^X, s < So, (18) 
where < g2{x^ Po) < 1 and g'^{x,PQ) > 0. 



^The primes denote diflferentiation with respect to the first argument. 
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Proof. We shall repeat the proof from [TT] in order to fix the notation for the rest 
of the paper. Relation so — s = giiro — i^, Po) implies 

d{sQ-s) dgi{ro~r,po) d{rQ ~ r) d{sa-s)/da 

—d^= diro-r) d{ro-r)/da=''^^^P'^- ^^^^ 

If 

d{so ~ s)/da 
0(^0 — rj/oa 

then one can easily see that 



/^i(«)=(^)' with r=^I^^l^-ll, for a>l. (20) 
From the first equation in (jlip . we have 



P / (a-l)(a7-i) a--l 
^ = V a + ^^=-- ■ (21) 



Therefore 



„// N 1 I ce a"' — 1 ^ ^ 1 



since a''' > 1 for a > 1 and 7 > 1. Using the fact that /'(a) > and (|2ip the 
Implicit Function Theorem yields that there exists a — ai{(3 / kpq) such that 

giiro-r,po)= /ii(a)U=ai(/3/Kp5) ■ (22) 

Jo 

Since 5i(/3,Po) = /ii(a), Si (/?, Po) = ' 3^ and ^ = Kp^-/'(a) > it remains 

dp da 

to prove that < hi{a) < 1 and < h[{a). From ([201) we have 



and 

0</.;(a) = 4.^|^.F', (23) 

since y > 1 and Y' > 0. The second part of the theorem can be proved using the 
same technique. □ 

Lemma 3.2. Let po < pi and (3/kp\ < 6 < /3/kpq. Then 

^ ^ ^ = ?I (24) 

We would need an estimate of the difference of Riemann invariants across two shock 
waves which is more precise than the one in It is provided by the following 
theorem. 

Theorem 3.3. Let < £ < ^, sq < Si, and take two Si curves originating at the 
points (tq^sq) = {po,uo) and (ro,Si) = {pi,Ui), which are continued to the points 
(r, s) and (r, 52)7 respectively. Then we have 

< (so - s) - (si - S2) < a x/i (ro - r) (si - sq) , (25) 

where C* is a constant independent of e, po and pi . 
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Proof. Let z'^ — sq — s, — si — S2 and w = ro — r (look at the diagram shown 
in Figure [1]) . 




(''o,'Sl) = (Pl>Wl) 



s 
k 



r 



Figure 1. Two 1-shock wave curves in r — s plane. 



By Theorem 13.11 and the Mean Value Theorem we know that for pi > po, there 
exists 9 such that 



dhi{a) 



da 



a=a{e) 



AO) (A -A] dfi^ (26) 

V i^pI kpI j 



where G ( ) . The definitions of hi and a imply 

'J^>0, ^>0 and A_A>o 
da do kPq Kpi 

so z° — > for pi > Po- We need to estimate the integrand in ((26|) . By (|20|l . we 
have 

dhi{a) da{e) ^ 4(y - 1)(7 + l)a^ 



da 



dO 



Thus, 



and 



i^Vl PoPi Jo 



«\/7 PaPi Jo 



Y - 1 
dp . 



-Me) 



d(3 (27) 
(28) 



From Lemma |3^ we know that da/dd > for (3/kp\ < 6 < /3/kpq. Hence, 



kpU 



<a{e) <a 



KPoJ 



(29) 
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Moreover, 



-^ = fia)< 2 J (^^^ = 2«^/^ for « . « f A 
Kpl V a — 1 \ Kp{ 

At this point, we use a majorization of z'^ — different from the one in [11] in 
order to obtain bounds for C* independent of e. By (f29|) and the above inequahty 
we obtain 




(30) 

^ U Corny* HoTO-'t', 



for a (t;^) — -'^ ^"^^ 7 > 1. Since da/dO > 0, it follows by 





^1 Pi 



Using 



p\- pI = — - So) 



and ^-^7 = together with (PT|) we finally get 



4(7 + 1) e 
z < — — (si - So) w 

4(7 + 1) ^ 



(32) 



%/7Po 

Suppose that p is the first component of the solution of the Riemann problem ([l] 
for uq > ui. Then Lemma 3.1 from [3] yields that for small e > 0, there exists 
C > independent of e, such that p < C/e. Now, using ((32)) . if po ~ 1/e then 
Po 1 as e ^ 0. For e small enough we may write Po > Ci. Thus, there exists a 
constant C* independent of k, po, pi and e such that 

- < a %/i (si - So) tw, (33) 

holds. This completes the proof of the theorem. □ 

The theorem that follows can be proved in the same way. 

Theorem 3.4. Let < e < tq > ti, and take two S2 curves originating at the 
points (ro,So) = {pq,uq) and (ri,so) = (pi,Ui), which are continued to the points 
{r,s) and (r2,s), respectively. Then we have 

< (ro - r) - (n - rs) < {so ~ s) (ro ~ ri) , (34) 
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where C** is a constant independent of e, po and pi. 

(Al) We shall use the following convention: C* denotes the maximum of the con- 
stants C* and C** from Theorems 13.31 and 13. 4( respectively. 

In the following theorem (3 and x denote 5*1 and 5*2, respectively, while o and tt 
denote Ri and R2, respectively. The prime is reserved for after interaction waves. 
(For example, the interaction of 5*2 and 5*1 which produces Si and ^2 is denoted 
hyx + P^P' + X'O 

Theorem 3.5. If < e < ^, then the following estimates are valid for the corre- 
sponding interactions: 

(1) ^2 and Si interaction: 

(a) x + P^f3' + x' 

l/3'l < Ix'l < \x\ + C.V^\mx\, or 

there exist 77, ^ such that 

(b) X + /? -> /5' + X' 

0< |/3'|-|/3|-e, Ix'l < Ixl + C.V^mxl+V, 
where < rj < g'i{\f3\, po)£. < or 

(c) p ^ f3' + ^' 

0< |x'l = Ixl-^, l/3'l < \f3\+aV^\xm+v, 
where < 77 < g'i{\x\,Po)C < C • 

(2) ^2 and Ri (or R2 and Si) interaction: 

(a) X + o ^ o' + x' 

Ix'l = 1x1, \o'\<\o\ + C.^e\x\\o\. 

(b) ^ + + V 

W\ = \(i\, Itt'I < H + av^i/^iM . 

(3) ^2 and 5*2 (or Si and Si) interaction: 

(a) Xi + X2 o' + x' : 

|x'l = IXi| + IX2|, |o'| < |xi| + |X2| • 

(b) Pi+(32^ (3' + n' : 

l/3'l = l/3i| + 1/321, k'l < |/3i| + |/32| . 

(4) S2 and R2 ( or Ri and Si ) interaction: 

(a) 1° X + ^ /3' + x' • there exist 1-shock (3q and 2-shock xo such that 
Ixol = Ixl - ^, \Po\=V and Xo + Po ^ + x' , 

where < 77 < SadxLpi)^ < ^ • 

2° X + TT — !■ + tt' .■ there exist rj, ^ such that 

k'|<M, l/3'l = 77 < C = Ixl, 

where < 77 < 52(lxl,Pi)C < £. ■ 

(b) 1° o + /3 — > /?' + x' : there exist 1-shock Po o,nd 2-shock xo such that 
\Po\ = \P\ - t Ixol = V and Xo + Po ^ P' + x' , 

where ()< < g'2{\Pl P2)£. < i ■ 
2° o + P^o' + X' 
|o'|<H, |x'|=77<e = |/3|, 
where < 77 < g'i{\Pl po)S. < C • 

(5) i?2 and S2 (or Si and Ri) interaction: 

(a) 1° TT + X ^ /3' + x' • there exist 77, ^ such that 
Ix'l = 1x1 \P'\^V, 
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where < ?/ < 5i(lxl,P2)^ < £. ■ 

W\<\n\, \f3'\ = < ^ ^ \x\, 
where < i] < gadxLPo)^ < C • 
(b) 1° /3 + o — !■ /3' + x' : there exist rj, ^ such that 

where < rj < g[{\l3\,pi)^ < £, . 
2° /? + o ^ o' + x' 
|o'|<H, \x'\^v<i^\P\, 
where < i] < g[{\l3\,pi)^ < £, . 
(6) i?2 and Ri interaction: 

TT + O ^ o' + Tt' 

W\ = \ol k'l = |^|. 
Here C* is a positive constant defined as in (Al). 

Proof. This theorem can be proved using the same tools as in [11^ and therefore 
wiU be omitted. The only differences are: the constant C* is now independent of 
s, (3,Xi POi Pi ^^^d p2, and we have instead of e in the estimates (l)(a) — (l)(c), 
(2) (a) and (2) (6). □ 

The main part of the paper is the interaction problem of delta shocks via pressure 
perturbation. Thus, one needs to control shock and rarefaction strengths as p goes 
to infinity as £ ^ (more precisely, when p is boimdcd by const/e). Because of 
that, we give their estimates in r — s plane based on Theorem 13.31 and Theorem 
13.41 Let {po,uo) = (fo,so) be connected with (p, u) = (''i s) by a 1-rarefaction 
(or 1-shock) wave, while (p, u) = ir,s) be connected with (pi,ui) = (ri,si) by a 
2-rarefaction (or 2-shock) wave. Then the strength of 1-rarefaction wave is 

2 

r -ro = ^{p^- p"), p<po, (35) 

and the strength of 2-rarefaction wave is 

2 

si - s = ^(pf -p^), p<pi. (36) 
The strength of 1-shock wave is estimated by 

P ^^o f P 



7/2 

2pgx/iln^ <ro-r<2-=^^ (il) -pg, p>po, (37) 
Po Vl + 2e \PoJ 

while, the strength of 2-shock wave is estimated by 

/- / \ 7/2 

2pfyiln-^<s-5i<2-=^ ^ -pl, p>pi. (38) 
Pi VI + 2e \pij 

Let us estimate the upper bound of the 1-shock wave given in (j37l) . For the function 
gi from (fTTl) we have < g'i{l3,po) < 1 and < gi{f3,po), so 

hm g'i(l/3|,Po)<l. 

Let us consider two special cases needed for our investigation. The first case: p > po 
and p ^ 1/e. We have that there exist constants C, C and C independent of e such 
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that 



2 -^y= • < 2 , ^ ■ -= < const. 

It follows that there exists a constant C2, independent of e and poi such that 

supgi(l/3|,Po) :-C2 < 1 . (39) 

H6I1C6 

The second case: p > pQ, p ^ 1/e and po ~ V^- Then 

^■p'^ const ^ 2 -=£= /^p- ^ 0(Vi) . 
Po Vl + 2e V /^o 

Again, — > cx3 is impossible and (UOl) holds. In order to estimate the strength of 
S2, we can use the same arguments to prove 

sup52(lxl,Po) C4 < 1, (41) 

and 

^^4^>^-C5>0. (42) 
92(\x\,Po) C4 

From now on, we shall put 

Co = min{C3,C5}. (43) 

4. Global interaction estimates 

This section contains all the necessary assertions from [1 with several changes in 
constants. All changes are similar to those from the previous section. 

Definition 4.1. [Tj A Lipschitz curve J defined by t ~ T(x), a; G R is called an 
I-curve, ij \T' {x)\ < 1/A. We denote J2 > Ji, if Ti 7^ T2 andT2{x) > Ti{x),x e K. 
Denoting by Sj{J) the set of j- shock waves crossing J and S{J) — Si{J) + S2{J), 
we define 

L-{J)^ J2 QiJ)- E l/^IW- (44) 

aeS(J) PeSi{J),xeS2iJ), I3,X approach 

Set F{J) = L^{J) + K ■ Q{J), where K := 4C*^/e. A space-like line lying between 
the initial line and the first interaction point is denoted with O. 

Lemma 4.2. 

Q(0)<i-(0)2. (45) 
Proof. The proof follows straightforward from Definition 14.11 □ 

Lemma 4.3. Assuming 4C* y/eL^{0) < I, we have 

F{0) < 2L-{0) . (46) 
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Proof. 

F{0) = L-{0) + k Q{0) < L-{0) + k L-{O f (by gS])) 

= L-{0){1 + kL-{0)) = L-{0){1 + 4 a VsL-{0)) 

< L-{0){1 + 1) = 2L-{0) . 



□ 



As in 2 , consider a interval C M and a map a : J ^ M". The total variation 
(TV) of a is then defined as 



TV{a) := sup { V \a{xj) ~ a(a;j_i)| 




where the supremum is taken over aU > 1 and aU (A^ + l)-tuples of points Xj G J 
such that xo < a;i < • • • < a; at. Now, we give a new estimate for L+(0). Here, 
i'^(O) denotes the sum of the rarefaction waves strengths which cross the hne O. 

Lemma 4.4. We have 

L-{0) < TV{ro{x),sa{x)) and L+{0) < TViMx), soix)) . (47) 

The estimates in previous Lemma can easily be verified. The uniform bounds of 
F{J) follows from the following theorem. 

Theorem 4.5. // V^FiO) < min |i, then ^(Jz) < F(Ji) for Ja > Ji. 

Particulary, (J) < F{0). 

Proof. This theorem can be proved in the same way as Lemma 5 from [llj and 
hence the proof will be omitted. One has just to substitute a constant K from the 
original proof with the determined value K here. □ 

Lemma 4.6. Assume that KL^{0) < 1 and that 

y^n/(ro(x),so(x))<i-.min|^,^Co| . (48) 

Then kF{0) < min{2, Co}. 

Proof. Using Lemma 14.31 and 14.41 we have 

^VeFiO) < V2^eL-{0) < TV {ro{x), so{x)) < ^.min | ^, ^ Coj 

Multiplying it with 8 C* / V2, one gets 

AC,V^F{0) = kF{0) < min{2,Co} . 
which proves the claim. □ 

The right hand side of does not depend on e, and then one can say that 
TV (ro(a;), so{x)) may be arbitrarily large since we can always choose e small enough 
in order to fulfill (|48|) with 7 = 1 + 2£. Then we can apply wave front tracking pro- 
cedure from [1] for each such e, and obtain a sequence of step functions converging 
to the entropic solution. One only needs to replace CeF{0) and from |lj with 
C*^F{0) and Co, respectively. 
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5. Approximate delta shock solutions to pressureless gas dynamics 

Our main task is to solve delta shock interaction problem for pressureless gas dy- 
namics model. Accordingly, we will introduce a solution concept from [3] (somewhat 
simplified) and check consistency of theoretical and numerical wave front tracking 
results by letting e — » 0. 



5.1. Basic notions. In this section we shall use the notions and assertions from 
[9]. It contains results for a 3 x 3 system with energy conservation law added, but 
all the results can also be applied to system ([2|), too. Let us start with the basic 
definitions. Vector valued function of the form 

'Uo, x<c{t)-ae{t) 
( )^,f/i.(0, c{t)-aAt)<x<c{t)_ 

^ U2,e{t), c{t)<X<c{t) + b,{t) ^ ^ 

Ui, X>c{t)+be{t) 

is called weighted shadow wave (weighted SDW, for short). Here, U := {p,u). The 
functions a^, are continuous functions satisfying a^{0) — xi^g and bg{0) — X2,s- 
The SDW is constant if Ui^e and [/2,e are just constants. If, in addition, xi.e ~ 
X2,s = 0, then the wave is called simple. 
The value 

aS) a,{t)Ui,e{t) + b,{t)U2At) 

is called the strength and c'{t) is called the speed of the shadow wave. We assume 
that lim^^o o'e(^) — '^{t) £ exists for every t>Q and 

lim / Ueix,t)(j){x,t)dxdt = (Uq + {Ui - Uo)9ix - c{t)) + a{t) 5{x - c{t)) , (t){x , t)) , 

for t > 0, where is a Heaviside function. The SDW central line is given by 
X — c(t), while X — c{t) — ae{t) and x — c{t) + be{t) are called the external SDW 
lines. The values xi^,, and X2,e are called the shifts, while Ui^e{t) and U2,e{t) are 
called the intermediate states of a given SDW. 

Let i e {1,2, ... ,n}. We assume ||?7*||loo = 0{e^^), if / and g have at most a 
linear growth with respect to i-th component, or otherwise ||C^e||L°° = o(e~^). The 
components of the first kind are called major ones, while the ones of the second 
kind are called minor ones. 

A delta shock is a SDW associated with a 6 distribution with all minor components 
having finite limits as e — *■ 0. 

The following lemma is the base of all calculations involving SDWs. 

Lemma 5.1. Let f,g E C{ft : M") and U : ft C K" be a piecewise constant 

function for every t > 0. Let us also suppose that f and g satisfy 

max{||/(;7,.e)||L», ||5(f/^,e)||L~} - 0{e-'). (50) 

1—1,2 
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Then 

{dtfiU,), 4>) « 1^ lim ^ (a,{t)f{UiAi)) + hS)f{U2At))) 0(^(0, t) dt 
c\t)(f{Ui)-f{Uo)) 4>{c{t),t)dt 



oo 



+ j limc'(t)(^a,(t)/(C/i,,(t)) + fe,(t)/(;72,eW)ja,0(c(i),i)rfi 

(51) 

and 

{d^g{U,),<j>) « / (g{U,) - giUo)) <P{c{t), t) dt 

^° 1 ^ ^ (52) 

lim ({aS)g{UiAt)) + {be{t)9{U2At))) dA{c{t),t) dt. 

5.2. Entropy conditions. Let ?y([/) be a semi-convex entropy function for ([2]), 
with entropy-flux function q{U). We shall use entropy condition in the following 
form. A weak or approximate solution [/^ = (p^, u^) to system ^ with initial data 
f7|j^Q = Uo^e is admissible provided that for every T > we have 

lime-.o / / ??(C/e)9t0 + 9(C/e)a,0dtda;+ /77(t^o,e(2:,O))0(:E,O)dx>O, (53) 
Jm Jo JR 

for all non-negative test functions (f> £ C^{M. x (— cx3,T)). 

Using Lemma 15.11 with / substituted by j] and g by g and the fact that the delta 
function is a non-negative distribution, the first condition for SDW Ug from 
to be admissible is given by 

-c'(i)(??(C/i) - viUo)) + {q{U,) - qiUo)) 

d (54) 

+ lim MU,At))ae + v{U2At))be) < 0. 
£^0 dt 

The derivative of delta function changes the sign, so [4 has to satisfy 
\imc'{t){r]{UiAt))ae+v{U2At))be) 

£— »0 

-q{UiAt)M) - q{U2At))be{t) - 

in addition. 

These conditions are much simpler in the case of simple SDW when Uq, Ui, Ui^e 
and U2.e are constants: 

Ih^,^o - civiUi) - 77(C/o)) + a^UiA + beViU2A + l{Ui) ~ q{Uo) < (56) 

and 

lim -c{aMUiA + beV(U2.e)) + a,q{UiA + b,q{U2A = 0- (57) 

In most of the papers with delta or singular shock solution, the authors use over- 
compressibility as the admissibility condition. A wave is called the overcompressive 
one if all characteristics from both sides of the SDW line run into a shock curve, 
i.e. 

X^{Uo)>c'{t)>X,{Ul), i = l,...,n, 
where c is a shock speed and x = \i{U)t, i = 1, . . . , n are the characteristics of the 
system. One will see that these notations coincide with our model case. 



(55) 
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The entropy condition is connected with the problem of uniqueness for a weak so- 
lution of the conservation law system. We give a definition of weak (distributional) 
uniqueness and some results about it afterward. 

Definition 5.2. An SDW solution is called weakly unique if its distributional image 
is unique. More precisely, a speed c of the wave has to be unique as well as the limit 

lim a^Ui e + 6e?72 e- 

Let i € {I, . . . ,n}. If a limit lim a^t/^ ^ + b^U2 ^ is unique, then we say that the i-th 
component is unique. 

Note that all minor components of are unique by default. 

5.3. Entropy solutions to Riemann problem for pressureless gas dynamics 
model. The proof for the following theorem in the case of 3 x 3 PGD model is 
given in [9J. Its restriction to a 2 x 2 system is straightforward and therefore not 
discussed here. 

Theorem 5.3. Suppose that uq > ui. Then there exists a unique shadow wave 
solution of the form J^ffp to the Riemann problem {H 0^ satisfying the entropy 
inequality i5S\) with rj and q as defined above. 

Moreover, the validity of \5S\) for all semi- convex entropies rj are equivalent to the 
overcompressibility of the shadow wave. 

Our aim is to show the structure of a solution in order to be able to compare it 
with a numerical approximation described above. For our purposes it is safe to take 
tte = be = e in the sequel. In the proof of Theorem 15.31 we showed that a SDW 
solution (with U — {p,u)) to and initial data ([3]), with uq > ui, had to 
satisfy 



if po 7^ pi, and c — Ug — (uo + ui)/2, if pQ — pi. That defines a weakly unique 
SDW solution to the problem. 

5.4. Two SDWs interaction. The main advantage of using weighted SDWs (in- 
termediate states vary with t in addition) is for solving SDW interaction problem. 
Then we can proceed with the main part of the paper by showing numerically that 
such a solution can be viewed as a limit of gas dynamics model with a vanishing 
pressure as perturbation. Note that verification of delta shock existence has already 
been obtained in [3j (see [Gj for a somewhat general model). 

Suppose that two SDWs interact in a point {X, T) . The superscript 1 is used for 
data in the left wave while the superscript 2 is used for the right one. The first 
SDW connects the states Uq = {po,uo) with Ui — {pi,ui), while the second one 
connects the states Ui — {pi,ui) with U2 — {p2,U2). 

Again, the following theorem has been proved in [9j for the extended PGD system, 
and the proof can easily be adopted for the present one ([2]). 

Theorem 5.4. The result of two SDW interactions for the pressureless system (0) 
is a weakly unique single entropic weighted SDW. 



c 



lim 




{us does not depend on e) 



lim epe 



c[p] - [pu] = {uq - ui)VPoPi, 
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We use the following notation: [x]i :— x\ — xo, \x\i :— xi — x\ and \x\ := X2 — xo- 
The weighted SDW solution from the above theorem satisfies the following: The 
speed is given by c'{t) = Us{t) :— lim^^^o Ue{t), while Us{t) and ^(t) := lime^o ^Peit) 
satisfies the following ODEs system 

{^{t)Us{t)y = Us{t)[pu] - [pu 

with the initial data 



(58) 



£,{T)u,(T) =(c C + c e )r = ^ji^ MiVPopT (5g<) 

N2\/pip2 r. 



[ph 

Here are some facts regarding the solution (^(t), Us(t)), t > T to the above initial 
data problem (see [S]): 

(1) ^(t), for t > r, is an increasing function when exists. The initial data 
£,{T) > and £,{t) is always positive function for t > T (when exists), since 

Uo> Ul> U2. 

(2) From the system ((58)) we have 

<{t) - -^iWlit) - 2[pu\uS) + W]). 

The value — l/^(t) is now always negative for t > T. The roots of the 
right-hand side of the above ODE are denoted as Ai < A2. Then, for 
[p] ^ 0, 

[pu] ± |uo - M2IVP0P2 

W ■ 

Assume that [p] > 0. If Us{t) e {Ai,A2), then Us(t) increases, and if 
Us{t) e (—00,^1) U (A2,+oo), then Ms(t) decreases. The opposite holds if 
[p] < 0. There are two possible cases: 

• If po > P2, then U2 < Ai < uq < A2. If Us{T) G (1*2,^1), then Us{t) 
increases for t > T but stays bellow Ai. If Us{T) e {Ai,uq), then Us{t) 
decreases for i > T but stays above Ai. 

• li p2 > po, then Ai < U2 < A2 < uq. Again, if Us{T) e {u2, A2), then 
Us{t) increases for t > T but stays bellow A2. If Us{T) G {A2,uq), 
then Us{t) decreases ior t > T but stays above A2. 

This implies uq > Us{t) > U2 (the SDW is overcompressive) . Also, one will 
see that numerical examples resemble these asymptotic properties of Us{t) 
as t ^ 00. 



6. Numerical results 

In this section one can find numerical results which show a consistency of theoretical 
(in the sense of SDWs) and numerical results. Consider system ([T]) with the initial 
data 

{{Pq,uo), X < ai 
(/9i,ui), ai < X < a2 (60) 
{P2,U2), X > a2 
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Table 1. Parameter description 



Parameter 


Description 


K 


Adiabatic constant defined in 


Pe 


First component of the intermediate state of the sohition for ([1] I60p . 


Ue 


Second component of the intermediate state of the solution for ([1] 16011. 


Cl 


Speed of the first left shock. 


C2 


Speed of the last right shock. 


\Eqi\ 


Left hand side of the integral on the first equation in ([1]). 


\Eq2\ 


Left hand side of the integral on the second equation in 



where ai < a2, wq > ui > U2- Then (see [3]), for e small enough, there exist 
(pi,£,iti,e) G M+ X M and (p2,e,W2,£) S M+ x K, so that: 

• (pojMq) is connected with (/Oi,e,Ui,e) by an 1-shock, and (pi,e,iti,6) is con- 
nected with by a 2-shock, 

• (pijMi) is connected with (p2.6,U2,e) by an 1-shock, while {p2,eiU2.e) is 
connected with (p2,M2) by a 2-shock. 

A numerical solution is obtained by wave front tracking algorithm described in [1] . 
In order to verify two delta shocks interaction, we shall consider two cases. 

Case A. Suppose that (poi^^o) is connected with (pi,Mi) by a single delta shock 
and (pijMi) is connected with {p2,U2) by a single delta shock, too. Assume that 
{POjUq) can be connected with {p2,U2) by a single delta shock (so-called simple 
SDW, see [9]). The resulting SDW has a constant speed as a consequence. That 
can be done by choosing a special value for p2 provided that po,uo, pi,ui and U2 
are already given. 

Case B. We choose arbitrarily p2, i.e. the resulting SDW has a variable speed (a 
central SDW curve is no longer a line). The numerical results are given in Tables 
HElandll 

6.1. Case A. 

Example 6.1. Let ai — 0, a2 — 2, {po,uo) — (1,1), (pi,ui) — (1.2,0.8) and 
U2 = 0.7. Now, for p2 = 1.14286, there exists a single simple SDW as a solution to 
the interaction problem. 

Table 2. (po, wo) = (1, 1), (pi, "i) = (1.2, 0.8), (p2, ^2) = (1.14286, 0.7) 



7 


K 


e 


Pe 


Ue 


Cl 


C2 


\Eqi\ 


\Eq2\ 


2 


0.5 


0.5 


1.29979 


0.80062 


0.13553 


1.53337 


2 • 10--^ 


3 • 10-^ 


1.2 


0.29 


0.1 


1.68612 


0.82806 


0.57746 


1.09746 


1 • IQ--^ 


2 • 10-^ 


1.02 


0.099 


0.01 


4.22718 


0.84163 


0.79256 


0.89411 


2 • 10-^ 


3 • 10-^ 


1.01 


0.071 


0.005 


6.71337 


0.84311 


0.81565 


0.87247 


8 • 10-^ 


2 • 10-^ 


1.006 


0.055 


0.003 


9.95729 


0.84379 


0.82636 


0.86254 


1 • 10-^ 


6 • 10-=* 



After interaction, the speed of the resulting wave is cs = 0.844994. Two SDWs will 
interact in a point {X, T) = (12.365, 13.8088) with such data. Now, we are going to 
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explain Figures [2| Ffl and 1121 which arc iUustrations of appropriate numerical results. 
For each s we have two piecewise linear half-lines. The left one originates from 
the point (x,t) = (ai,0), while the right one originates from the point (x,t) = 
(02,0). The i-th linear segment of these half- lines can be written in the form 
X = Cij (t — ti) + Xi, i > 1, j = 1, 2, Xi < X < Xi+i, ti < t < ii+i, where Ci^i stands 
for the speed of the first (Si) wave on the left hand side in phase plane, while Ci^2 
stands for the speed of the last {S2) wave on the left hand side in phase plane at 
each i-th segment. Interactions of the waves occur at the points {xi,ti), i > I. 
After two delta shock interaction, the resulting delta shock central line in Figure [2] 
(dashed fine) starts from {X,T) and it is calculated exphcitly from system ([2]). 

6.2. Case B. 

Example 6.2. Let ai — 0, a2 — 2, {po,uo) — (1,1), (pi,ui) = (1.2,0.8) and 
(P2,W2) = (1.3,0.7). Two SDWs will interact in a point {X,T) = (12.2291, 13.657) 
with such data. After two delta shock interaction, the resulting delta shock central 
lines in Figures [7l and [T2l f dashed lines) start from {X,T), too. Here, 

t 

x{t) = J Us{p) dp + X, t>T, 

T 

while Us{t) represents the second component of the solution {£,{t), Us{t)) of system 
([55)1 with initial conditions 



Table 3. (po, uo) = (1, 1), (pi, wi) = (1-2, 0.8), (p2, ^2) = (1.3, 0.7) 



7 


K 


e 


Pe 




Cl 


C2 


\Eqi\ 


\Eq2\ 


2 


0.5 


0.5 


1.38131 


0.74968 


0.09312 


1.54395 


3 • 10-* 


5 • 10-4 


1.2 


0.29 


0.1 


1.79287 


0.80661 


0.56268 


1.08778 


5 • 10-* 


6 • 10-* 


1.02 


0.099 


0.01 


4.49667 


0.83356 


0.78596 


0.88787 


3 • 10-^ 


3 • 10-^ 


1.01 


0.071 


0.005 


7.14038 


0.83646 


0.80983 


0.86684 


1 • 10-^ 


3 • 10-^ 


1.006 


0.055 


0.003 


10.5884 


0.83782 


0.82091 


0.85711 


3 • 10-^ 


2 • 10-^ 



Example 6.3. Let ai — 0, a2 — 2, {po,uo) = (1,1), (pi,ui) ~ (0.8,0.9) and 
(/02, U2) = (0.9, 0.7). Two SDWs will interact in a point {X, T) = (12.2364, 12.8427) 
with such data. 



Table 4. (po, uo) = (1, 1), (pi, ^i) = (0.8, 0.9), (p2, ^2) = (0.9, 0.7) 



7 


K 


e 


Pe 




Cl 


C2 


1^91 1 


\Eq2\ 


2 


0.5 


0.5 


1.16621 


0.88674 


0.20529 


1.51813 


8 • 10--^ 


4-10-^ 


1.2 


0.29 


0.1 


1.50419 


0.86711 


0.60356 


1.11606 


2 • 10--* 


1 • lO-'i 


1.02 


0.099 


0.01 


3.75748 


0.85659 


0.80459 


0.90592 


4 • 10-'^ 


2 • 10-=* 


1.01 


0.071 


0.005 


5.96447 


0.85544 


0.82632 


0.88306 


1 • 10-'^ 


8 • 10--* 


1.006 


0.055 


0.003 


8.66370 


0.85341 


0.83808 


0.86341 


2 • 10-^ 


3 • 10-^ 
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7. Appendix 
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Figure 2. Phase x — t plane, Case A, Example 16.1 
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Figure 3. Speed of delta shock formed after double delta shock 
interaction, Case A, Example 16.11 
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Figure 4. Speed of the first {Si) wave on the left hand side and 
the last (52) wave on the right hand side for various e, Case A, 
Example 16.11 
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Figure 5. Solution p{x,t) for various e at t = 15000, Case A, 
Example 16.11 
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Figure 6. Solution u{x,t) for various e at t = 15000, Case A, 
Example 16.11 
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Figure 7. Phase x — t plane, Case B, Example 16.21 
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Figure 8. Speed of delta shock formed after double delta shock 
interaction, Case B, Example 16.21 
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Figure 9. Speed of the first {Si) wave on the left hand side and 
the last (5*2) wave on the right hand side for various e, Case B, 
Example 16.21 
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Figure 10. Solution p{x,t) for various e at t = 15000, Case B, 
Example 16.21 
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Figure 11. Solution u{x,t) for various e at i = 15000, Case B, 
Example 16.21 




Figure 13. Speed of delta shock formed after double delta shock 
interaction, Case B, Example 16.31 
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Figure 14. Speed of the first (5*1) wave on the left hand side and 
the last (5*2) wave on the right hand side for various e, Case B, 
Example 16.31 
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Figure 15. Solution p{x,t) for various e at t = 15000, Case B, 
Example 16.31 
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Figure 16. Solution u{x,t) for various e dX t — 15000, Case B, 
Example 16.31 



